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We have investigated two-electron bound state formation in a square two-dimensional t — J — U 
model with hopping anisotropies for zero electron density; these anisotropies are introduced to mimic 
C ' the hopping energies similar to those expected in stripe-like arrangements of holes and spins found 

in various transition metal oxides. In this report we provide analytical solutions to this problem, 
and thus demonstrate that bound-state formation occurs at a critical exchange coupling, J c , that 
decreases to zero in the limit of extreme hopping anisotropy t y /t x — > 0. This result should be 
5^ , contrasted with J c /t = 2 for either a one-dimensional chain, or a two-dimensional plane with 

■ isotropic hopping. Most importantly, this behaviour is found to be qualitatively similar to that of 

two electrons on the two-leg ladder problem in the limit of tmtcrchamAmtrachain — * 0. Using the 
latter result as guidance, we have evaluated the pair correlation function, thus determining that the 
bound state corresponds to one electron moving along one chain, with the second electron moving 
along the opposite chain, similar to two electrons confined to move along parallel, neighbouring, 
metallic stripes. We emphasize that the above results are not restricted to the zero density limit — 
we have completed an exact diagonalization study of two holes in a 12 x 2 two-leg ladder described 
by the t — J model and have found that the above-mentioned lowering of the binding energy with 
hopping anisotropy persists near half filling. 
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\^ • The experimental verification for the existence of stripes in layered transition metal oxides, such as the high-T c 
superconductors, is now robust Q, and yet an important question remains: Does the presence of such rivers of charge 
help, hinder, or even possibly create, the pairing instability that leads to superconductivity? Simple model calculations 
that might shed light on this question are clearly of value, and in this paper we present one such study for a model 
Hamiltonian that we have considered in an attempt to mimic some of the physics of the stripe phases. 

Previous work on the magnetic properties of the very weakly doped cuprates |2j modelled the observed experimental 
support for stripe correlations jj| using an effective Hamiltonian in which a (spatially) anisotropic exchange interaction 
was implemented to represent the stripe-induced magnetic energy scales. That is, in the direction parallel to the stripes 
the full local Cu — Cu exchange would be present, while perpendicular to the stripes a reduced exchange would be 
I ■ encountered across such rivers of charge. Renormalized Hamiltonians of a similar simplifying spirit were also used in 
other studies of the doped cuprates ^|(| . 

To this end, we have considered an effective Hamiltonian in which a carrier's ability to move along the river of 
charge is far larger than its freedom to move between the rivers. This reduced mobility reflects a carrier's reduced 
ability to move in an antiferromagnetically correlated region, such as exists between the rivers of charge. The simplest 
theoretical modelling of such a situation corresponds to a Hamiltonian with local magnetic exchange that is the same 
k> , in all directions, but with a hopping parameter which is very different in one direction than the other. (Several studies 
• of holes moving on 2-leg ladders near half filling have considered a similar problem, but all but one have employed a 
& \ limit that is different to that which we have investigated. To be particular, these authors have considered the hopping 
between the chains of the ladder to be larger than that along each chain 0-0. The recent work of Weihong, et al, 
jllj considers a limit similar to that of this paper, although their focus is two holes near half filling, not two electrons 
near zero filling.) 

As a consequence of the above considerations, in this paper we consider the so-called t — J — U model, which is 
defined by 

H = - ^ M c l,a c ^ + h.C.) + ^JijfSi • Sj - 2 n i n j) + U Y1 ni <t n U ' W 
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where the sites of a two-dimensional square lattice of size L x x L y with periodic boundary conditions are labeled by 
the indices i and j, t^j and are the hopping integrals and exchange couplings between sites i and j, respectively, 
Cj lCT is the annihilation operator for an electron of spin a, U is the on-site Hubbard repulsion energy, and nj )(T is the 
number operator for electrons at site i with spin a. 

The most familiar strong-coupling variant of the Hubbard model is the t—J model, and the physics of (square lattice) 
doped Mott insulators described by this model was reviewed by Dagotto JT^l . As emphasized by, e.g., Anderson [fl3|| , 
the vital component of the t—J Hamiltonian is the constraint of no double occupancy. That is, in the t — J model 
one does not use the electron creation and annihilation operators of Eq. ([!]) , but rather one uses constrained creation 
and annihilation operators (for example, see the discussions in Ref. [Q). However, the above t — J — U Hamiltonian 
can be used to accomplish this same mathematical projection by taking U — > oo — this simplifying approach has 
been noted by a variety of researchers (see, e.g., Refs. |i~4|-|l6[] ) , and will be used by us also. 

In this paper we restrict t%j and Jy to be nonzero for near neighbours (NN) only. Further, we allow the hopping 
integral in the x direction, t x , to be different than the hopping integral in the y direction, t y . We have also investigated 
the physics that arises when J x is allowed to be different than J y , but find that no qualitatively new physics arises 
as long as both J x and J y remain nonzero. Thus, from now on we set J = J x = J y , and we analyze the resulting 
hopping anisotropy problem in terms of 

t x = t , r= t -f . (2) 

By taking the limit U/t — > oo the only two (dimensionless) energy scales left in this problem are J/t and r. 

The physics of the t — J —U model for isotropic NN hopping (r = 1) in the dilute electron density limit is clear: the 
antiferromagnetic exchange provides an attractive interaction between the electrons on neighbouring sites forming a 
singlet state; however, a bound state that is formed from this interaction must do so in a manner that excludes double 
occupancy, a result that must follow in the U — > oo limit. As mentioned above, this problem for isotropic hopping 
and exchange has already been studied elsewhere JlJJl^]. Here we shall analyze the bound state exchange threshold, 
to be denoted by J c /t, as a function of hopping anisotropy r. 

We consider the situation in which there are only two electrons present, in which case we can write for the singlet 
wavefunction of a two-electron system as ]14|Jl7| 

ip = J2^j)4Ai\°) ( 3 ) 

where |0) denotes the empty lattice. Any bound state can be solved for by the direct solution of the Schrodinger 
equation. Expressing 4>(i,j) in terms of its Fourier transform, the latter of which involves the centre of mass and 
relative momenta, given by Q = ki + k2 and q = §(ki — ka), with <f>(ki, k2) = 4>q{<1), where 

0(k 1; k 2 ) = ±Y / e- ik *-''- ik *-'>4>(i,j) (4) 
with N = L x L y , the resulting equation for the bound-state energy E and its associated wave function 0Q(q) is 

, ( n) ^E k Q (k)-^E k J(q-k)0 Q (k) 

•Ml) = ~ tq- — ; 7q : • ( 5 ) 

E-e[f +q) -e(f - q) 

The dispersions e(k) and J(k), the former of which includes the hopping anisotropy, are given by e(k) = — 2t(cosk x + 
r cos k v ) and J(k) = 2J(cosfc2; +co$k y ) . 

The lowest energy bound states for the t — J — U model has zero centre-of-mass momentum (Q = 0). Then, one 
can solve Eq. (g) via the 3 x 3 set of algebraic equations: 

Co = UIqCq — 2JI X C X — 2JI y C y (6) 
Cx = UI x Cq 2JI xx C x 2JI xy C y 

Cy = U IyC() 2JI xy C X 2JI yy Cy 



where 



Co = ^E^ k ) G^-^cobMW C, = l^cosfc^(k) (7) 

k k k 
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and 



E + 4i(cos q x + r cos q y ) 



COS< 



q 

Iy = -Y 

y N 

/ 

11 AT ^ £7 + 4<(cosg x + rcosq y ) 



COS Qj; 



U + 4£ (cos g x + r cos q y ) 



q 



E + 4i(cos ? x + r cos g y ) 
cos q x cos q a 



(8) 



AT ^ E + 4i(cos q x + r cos q y ) 



COS 2 q x 



lyy- N 



/v 



cos z g a 



-E + 4<(cos + r cos q y ) 



A methodology leading to the solution of such a set of equations is presented in the papers of Lin and Pethukov, 
et al. pi However, in those papers the authors were considering a two-dimensional square plane with isotropic hopping, 
and thus were able to exploit C x = C y owing to the result that the ground state wave function has s-wave symmetry. 
In our problem the Hamiltonian has a reduced symmetry, leading to C x ^ C y (for r ^ _1) , and thus we have employed 
a different approach — we will explain this method in detail in a future publication Jl7|] . 

We find that the resulting bound state is always symmetric under 7r/2 rotations, and reflections along the x or y 
directions; that is, the ground state transforms under the identity representation of the rectangular point group, and 
thus continuously interpolates with the s-wave bound state for r — 1. The threshold J c /t as a function of r is given 
by the following expression: 



t 



where / and g are functions given by 

4(r 3 - r 

f(r) 



167r 2 r 5/2 

g( r ) 



■ r — 1) sin 



/(r)-\//(r) 2 - 



-2nr 3 + 4r 5 / 2 + 2nr 2 + 4^? 
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(10) 



j(r) = 2(r 2 - 1) sin" 1 - ^ 2 + 2r 3/2 + 2V^ 



The resulting J c /t vs. r is plotted in Fig. Jlh from this plot we note (i) the r = 1 limit of J c /t = 2 result is found, 
confirming the known isotropic result [ |14}|lq |, and (ii) the surprising result that as r — > 0, J c /t — > 0. This latter 
result is important — it implies that one can form bound states for arbitrarily small J/t if one is able to create a 
situation in which the hopping anisotropy is large enough. Further, it implies that for a fixed J/t, the binding energy 
of two electrons will become increasingly larger as this anisotropy is increased. We have plotted this bound state 
energy result in Fig. g, for J/t — 2, from which it is seen that one can obtain remarkably large binding energies. 
While acknowledging our inability to relate this system directly to the high-T c systems, without further trepidation 
we evaluated the r — > limit of the binding energy for J/t = 0.3. We found a binding energy of 0.0068i, and in units 
in which t is roughly an eV, this implies a binding energy of roughly 100 K, a provocative energy. 
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FIG. 1. The dashed line is a plot of Eq. (|9j), and shows how the threshold exchange for bound state formation, J c /t, varies 
as a function of hopping anisotropy, r, for an infinite two-dimensional plane. The solid line is a plot of Eq. (||), which is the 
same quantity for a two-leg ladder. Note the rapid decrease to J c /t = for both systems as r — > 0. 
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FIG. 2. The value of the bound state energy, Ebs, relative to the bottom of the two-electron noninteracting band, which 
is equal to — 4t(l + r), as a function of r for J/t = 2 for the infinite two-dimensional plane. (Note that for this value of J/t a 
bound state is found for all r.) 



The question remains, what are the nature of the correlations in such bound states? In an attempt to answer this 
question we have solved the same problem as above not for an infinite two-dimensional plane, but for a two-leg ladder. 
Technically, this involves nothing more than using 

5>*E if]*** (id 

q ? B =0,7r J 

in the wave vector sums that appear in the evaluation of the various (elliptic) integrals of Eqs. (@,|[). The resulting 
J c /t as a function of r is given by 

j ( (1 + r) 2 + 2V^T^ - V(l + r 2 ) 2 + Ar) 

T = " TT 1 ( 12 ) 

t r + 1 

and is also shown in Fig. [j]. Thus, we see that both the infinite two-dimensional plane and the infinitely long two-leg 
ladder behave in a similar fashion. Further, the non analytic behaviour of J c /t is quantitatively similar in the extreme 
hopping anisotropy limit: in the limit of r — > the two-dimensional infinite plane has J c /t ~ (37r/4)y / r, while the 
two- leg ladder has J c /t ~ 2y / r. 

We have exploited this similarity by evaluating the pair correlation function for the two-leg ladder, noting that the 
simpler geometry of the two-leg ladder more readily allows for us to understand the spatial character of the extreme 
hopping anisotropy bound states. More specifically, we have calculated the probability that the two electrons are 
on the same chain or are on opposite chains. Our results as a function of r are shown in Fig. |3|, from which it is 
seen that the character of the small r bound states is two electrons moving on opposite chains, with the pairing 
resulting from a spin exchange interaction between the chains (when the electrons are on neighbouring sites). There 
is a higher energy bound state that corresponds to J c /t = 2, but in this bound state the systems behaves like two 
one-dimensional chains, and the threshold for this higher energy bound state is the same as for a one-dimensional 
t — J — U(U — > oo) model, viz. J c /t = 2. Thus, the lower energy and lower threshold bound states can only be 
achieved when the electrons are moving along opposite chains. 
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FIG. 3. The pair correlation for the bound state of the two-leg ladder was evaluated for all electron positions, and then 
summed to yield the probabilities that the electrons are on the same chain or on opposite chains. This quantity is shown as a 
function of hopping anisotropy, r. Note that as r —* the electrons are restricted to be on opposite chains only. 



This situation is qualitatively similar to two electrons moving in neighbouring stripes that still experience an 
exchange coupling between the stripes due to the antiferromagnetically correlated domain between the stripes, but 
which have an effective hopping frequency that is much larger along the metallic rivers of charge than between such 
rivers. In this very simplistic picture, it seems that stripe correlations leading to parallel rivers of charge separated 
by antiferromagnetically correlated domains can (at least) enhance bound state formation. 

We have been able to complete the above work analytically because of the relative simplicity of the zero density 
two-electron problem. Of course, the natural question is: Does this extrapolate to systems that are doped away from 
half- filling? To address this question we have considered the problem of two holes in an otherwise half-filled 12 x 2 
two-leg ladder for J/t = 0.3. For all r the ground state is Q = and transforms under the identity representation 
of the rectangular point group. The binding energy is again found to be substantially enhanced by the hopping 
anisotropy, a result previously noted in Ref. |ll| — our data is listed below in Table I. Lastly, we again evaluated the 
probability that the holes existed on the same or opposite chains of the two-leg ladder, and our results are qualitatively 
similar to those of the two-electron problem that we have shown in Fig. [| (the two-hole correlations near r = 1 have 
a greater probability of being on the opposite chains (about 0.62), but otherwise these two systems behave in very 
similar ways). 

TABLE I. The two hole binding energy, Eb, defined as Eb = E^h — "2E\h + Eoh, calculated for a 12 x 2 two-leg t — J ladder 
with J/t — 0.3, as a function of anisotropy parameter r. 



_r E b 

1.0 -0.1189 

0.8 -0.1614 

0.6 -0.2224 

0.4 -0.2803 

0.2 -0.3137 

0.1 -0.3254 
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Thus, the zero density two-electron problem and the half filling two- hole problem behave in a very similar manner, 
and are suggestive of the benefits of stripes in forming bound states. Whether or not this physics extrapolates to 
pairing instabilities and, e.g., enhances the mean-field superconducting transition temperatures, is being investigated 
presently. 
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